X-ray spectrum of a pinned charge density wave 



O 

o 

(N 



Alberto Rosso and Thierry Giamarchi 
Universite de Geneve, DPMC, 24 Quai Ernest Ansermet, CH-1211 Geneve 4, Switzerland 

We calculate the X-ray diffraction spectrum produced by a pinned charge density wave (CDW). 
The signature of the presence of a CDW consists of two satellite peaks, asymmetric as a consequence 
of disorder. The shape and the intensity of these peaks are determined in the case of a collective 
weak pinning using the variational method. We predict divergent asymmetric peaks, revealing 
the presence of a Bragg glass phase. We deal also with the long range Coulomb interactions, 
concluding that both peak divergence and anisotropy are enhanced. Finally we discuss how to 
detect experimentally the Bragg glass phase in the view of the role played by the finite resolution 
of measurements. 
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PACS numbers: 

The study of disordered elastic objects sheds light on 
the physics of a wide range of systems. A first class of 
systems is accurately modeled by elastic manifolds in the 
presence of randomness; significant examples include do- 
main walls in magnetioi or ferroelectrio^ materials, con- 
tact lines of liquid menisci on a rough substrate^- and 
propagating cracks in solids'*. A second class of dis- 
ordered elastic systems is given by periodic structures 
such as charge density waves (CDWs)^, vortex flux lines 
in type-II superconductors^ and Wigner crystals!. It 
was recently shown that periodic systems have unique 
properties, quite different from the ones of the inter- 
faces. In fact, if topological defects (i.e. dislocations 
etc.) in the crystal are excluded, displacements can grow 
only logarithmicall}iSi2ii£, in contrast with the power- 
law growth of interfaces. The positional order is only 
algebraically destroyediSiii, leading to divergent Bragg 
peaks and a nearly perfect crystal state. Quite remark- 
ably, it was shown that for weak disorder this solu- 
tion is stable to the proliferation of topological defects, 
and thus that a thermodynamically stable phase hav- 
ing both glassy properties and quasi-long range posi- 
tional order exists^i. This phase, nicknamed Bragg glass, 
has prompted many further analytical and experimen- 
tal studies (see e.g. Ref. [l^l3| for reviews and further 
references). Although its existence can be tested indi- 
rectly by the consequences on the phase diagram of vor- 
tex flux lines, the most direct proof is to measure the 
predicted algebraic decay of the positional order. Such a 
measurement can be done by means of diffraction exper- 
iments, using either neutrons or X-rays. Neutron diffrac- 
tion experiments have recently provided unambiguous 
evidencesi^ of the existence of the Bragg glass phase for 
vortex lattices. 

Another periodic system in which one can expect a 
Bragg glass to occur are CDWs^, where the electronic 
density shows a sinusoidal modulation. As a conse- 
quence of the electron-phonon interaction, this modula- 
tions generates a permanent distortion of the underlying 
lattice. This distortion can be revealed thanks to X-ray 
measurements: in fact, the corresponding X-ray spec- 
trum presents satellite peaks around each principal Bragg 
peak. These satellites contain information concerning the 



positional order of the CDW. In particular, we are inter- 
ested in the detection of effects due to disorder—. The 
X-ray experimental resolution is in principle much higher 
than the one that can be achieved by neutrons for vor- 
tex lattices, consequently CDW systems should be prime 
candidates to check for the existence of a Bragg glass 
state^^. However, compared to the case of vortex lattices 
the interpretation of the spectrum is much more compli- 
cated for two main reasons: (i) the phase of the CDW 
is the object described by an elastic energy, whereas the 
X-rays probe the displacements of the atoms in the crys- 
tal lattice (essentially a cosine of the phase); (ii) since 
the impurities substitute for some atoms of the crystal, 
the very presence of the impurities changes the X-ray 
spectrum. This gives rise to non-trivial terms of interfer- 
ence between disorder and atomic displacementsiSiii. It 
is thus necessary to make a detailed theoretical analysis 
of the diffraction due to a pinned CDW. 

In the past, the study of the spectrum has been 
carried out only either for strong pinning or at high 
temperaturesiSiiLiSii^. In this paper we focus on the 
low temperature limit, where a well formed CDW exists, 
and on weak disorder, for which one expects to be in 
the Bragg glass regime. Both the short and long range 
screening of the Coulomb interactions are considered. We 
show that in both cases the diffraction spectrum consists 
in two asymmetric peaks. The peaks are power-law di- 
vergent, with a stronger anisotropic shape in the case of 
unscreened long-range Coulomb interaction. This find- 
ing is consistent with the Bragg glass behavior—. The 
asymmetry divergence follows a subdominant power-law 
as well, with an exponent that we determine. A short 
account of part of the results of this paper was published 
in Ref. 

In section I we derive the model used to describe the 
interaction between the CDW and impurities. Two elas- 
tic limits are considered: if free electrons are present the 
elasticity has a simple short-range form, while, in the un- 
screened case. Coulomb interactions are responsible for 
a long range strongly anisotropic elastic term. In sec- 
tion II we discuss the X-ray intensity spectrum behavior 
in presence of a pinned CDW. In particular, we derive 
the different contributions to the satellite peaks and we 
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study their symmetry properties. In section III we evalu- 
ate explicitly the different terms by means of the replica 
techniques. Section IV contains the physical discussion 
in view of all the results obtained in this paper. The 
reader not interested in the details of calculations may 
move directly to this section skipping the previous one. 
Finally, in appendix A we evaluate the triplet contribu- 
tion and in appendix B we calculate the function [a] , used 
in section IV. 



I. THE MODEL 

The system we have studied is a CDW in a three di- 
mensional space. The electron density has the form^^: 



Pir) = Po + 



-V0(r) + pi|V|cos(gr + (^(r))), (1) 



where a single sinusoidal deformation of modulation vec- 
tor Q is considered and tp — \ip\e'^'^ is the CDW order 
parameter normalized to unity at T = 0. The first term 
of Eq. (Q), po, is the average density. The second one cor- 
responds to a density averaged at scales larger than Q~^. 
This contribution, also called forward scattering, encom- 
passes the local changes of the electron density related to 
the compression modes. The last term, also called back- 
ward scattering, describes the sinusoidal modulation at 
a scale of the order of Q~^. We can neglect all contri- 
butions stemming from higher harmonic terms as they 
are known to be important only at very low tempera- 
ture. The effective Hamiltonian can be obtained by a 
Ginsburg-Landau expansion of the order parameter 



am' 



(2) 



where a = [Tc — T) /Tc and 6 is a parameter whose value 
is defined by the microscopic theory. The configuration 
at minimum energy corresponds to — ^/a with the 
phase 4> equal to a constant. Around this equilibrium so- 
lution, fluctuations involve both the amplitude and the 
phase of the order parameter. We remark that while the 
first ones are more expensive in terms of energy varia- 
tions, due to the presence of the quadratic term, the sec- 
ond ones are massless. Following the model developed by 
Fukuyama-Lee-Rice (FLR)MiS24S, we take into account 
only the phase fluctuations (i.e. we neglect amplitude 
fluctuations). The mass of these excitations, called pha- 
sons, turns out to be quite large, because it depends on 
the ionic mass via the electron-phonon interaction. For 
this reason, we can neglect in the Hamiltonian the ki- 
netic term giving rise to quantum fluctuations. Within 
these hypotheses, the elastic Hamiltonian associated to 
the CDW reads 



(3) 



where c = U^ij}]"^ is the elastic constant. The thermal 



fluctuations of the phase, denoted as 
uated: 



are easily eval- 



2T 



\x)) = J V<j>e 



(4) 



The integral in q extends all over the Brillouin zone. In 
the low temperature regime, the fluctuations are thus 
small enough to guarantee the presence of an ordered 
phase. Moreover, we observe that the value of the in- 
tegral Q is related to the shortest length scale in the 
problem. Here we assume that this cut-off momentum, 
A, is given by the periodicity of the CDW (A Q)^. 

The form of the Hamiltonian |j2Jl is actually anisotropic 
along the Q-direction. In fact, a compression along Q 
produces an increase of the electric charge density which 
yields an increase of the stiffness, whereas all distortions 
along the other two directions do not involve any change 
in electrostatic energy. 

We evaluate the contribution of Coulomb interactions 
screened beyond the characteristic length A. Without 
any loss of generality we assume Q || a; and y is along one 
of the other two equivalent directions. The electrostatic 
energy takes the form: 



C/cx 



1 

2V 



d^rd^r'e 



-\r-r'\/\ Pir)p{r') 



(5) 



The main contribution to the electrostatic energy comes 
from the variations of the electron density p, expressed 
by Eq. Therefore, we restrict to consider only the 
forward scattering term: 



U (X 



pM' 
pI\M' 



d^,d^,',-\r-rvMr)d.^{r') 
\r — r'\ 

[qA? 



(6) 



This term introduces a g-dispersion in the elastic con- 
stant 



Cliq) cx 



BZ 

A2 



\m\' 



1 + g2A2 



Two regimes can be identified as a function of the screen- 
ing length A. (i) In the first one, valid for A not very large, 
we can neglect the dispersion in q and the resulting elas- 
ticity is short-range. The effect of Coulomb interaction 
is an enhancement of the elastic constant along the x- 
direction 



i/ei = / dxd'y ^ (9,0) 



Cl 



(7) 



By redefining the spatial variables. 



= x/^/ci and 



y/\/c, with c = (cic2)2j the Hamiltonian jT)) can 
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FIG. 1: Dispersion relation of the phase mode given by 
The transversal mode {9 = 90°) is acoustic, the longitudinal 
mode is gapped = 0°). The gap is reduced by the presence 
of longitudinal components (6 — 60°). 



be finally turned into an isotropic form, (ii) The sec- 
ond regime is characterized by large values of A. In 
this case, the electrostatic energy takes the form U ~ 
lBZ^^x/l^)\^il)\^ ^ ^iid consequently a long-range term 
appears in the elasticity: 



Hci. — 



dqx (Pq 
27r (27r)2 



2 2^ 



i0(9)r 



(8) 



At this stage, we briefly discuss the dispersion relation 
of different elastic regimes. The full FLR Hamiltonian is: 



FLR 



BZ 



-^PP 

2M " ' 



2 



(9) 



where Pq is the Fourier transform of the momentum den- 
sity and M is the phason mass density. The first term 
gives the kinetic energy and the second the elastic en- 
ergy. Using the standard canonical transformation we 
derive the corresponding dispersion relation: 



where 5](r) is the impurity probability distribution. Long 
range interactions are neglected and Vq is a positive con- 
stant which measures the impurity potential. At last the 
sign -I- (— ) is related to the repulsive (attractive) interac- 
tion between the electrons and the local impurity. Above 
two dimensions we can drop out the forward scattering 
term in the development of p{r). In fact, this term leads 
only to a trivial redefinition of the correlation functionsii. 
If Vo is small (the opposite case, the effect of a strong 
impurity, is discussed in Ref. [23) the FLR model for 
the elasticity is justified. In this case, the collective pin- 
ning S(r) is well described by a Gaussian distribution 
with zero average (we can always incorporate the effect 
of the averaged disorder into the bare parameters) and 
the correlator is given by I](r)E(r') = NjS{r — r'), with 
Nj — nj(l — nj). In the following, we will restrict our 
analysis to the repulsive case: pi\4'\ is absorbed in Vq and 
we define the disorder strength D = Vq'^Nj. 

On one hand, the disorder favors local distortions of 
the phase cj), on the other hand these deformations in- 
crease the elastic energy. A natural size Ra is defined if 
we consider the region where (fi varies by 27r. A simple 
energetic balance gives for a d-dimcnsional CDW 



„ 27r 

Etnt — - 



2 \ Ra 



(12) 



Optimizing the gain in potential energy versus the cost 
in elastic energy we get Ra — {c^ / D)^/'^^''^^ (for d — ?, 
Ra — c^/D). This length, called Fukuyama-Lee length 
(or Larkin-Ovchinikov length)iS*2&, is interpreted as the 
correlation length of the system. In this scenario the 
equilibrium state is always disordered and the long-range 
coherence in the phase is lost. Nevertheless, in this paper 
we show that, as in the case of vortex lattice, the latter 
prediction is correct only if we consider scales smaller 
than Ra, while it breaks down at larger distances. In 
particular, we see that it is the Fukuyama-Lee length to 
define the crossover between the short distance regime 
and the asymptotic one. 



uj{q) 



2isl 

M 



(10) 



If we decompose the vector q in its longitudinal {qx — 
qcosO) and transversal {q± = qcoaO) components, it is 
clear that in the short-range case q obeys a linear disper- 
sion law, with a slope equal to a/ (ci cos(6') -|- csm{6))/M. 
The dispersion for the long-range elasticity (jS)) is dis- 
played in Fig.^ the transversal modes remains acoustic, 
the longitudinal ones instead develop a gap^**. 

Finally, we consider the effect of a distribution of im- 
purities with concentration nj. The simplest coupling 
with the electron density is expressed by: 



Hdis = ±1^0 / rfVE(r)p(r) 



(11) 



II. SPECTRUM INTENSITY 

The X-ray diffraction is a powerful tool to detect 
any subtle change of the perfect crystalline structure. 
The electron density modulation is accompanied, via 
the electron-phonon interaction, by a lattice distortion 
u given by 



u{r) 



Uq 



V cos(Qr -t- (j){r)) cx W p{r). 



(13) 



Thus, the CDW instability produces a permanent sinu- 
soidal displacement of the atoms from their equilibrium 
position. This deformation is signaled, in the X-ray spec- 
trum, by the presence of satellite peaks around each prin- 
cipal Bragg peak. The analysis of the shape, the intensity 
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the disorder and over the thermal fluctuations, fi rep- 
resents the total amplitude scattered by the atom at the 
position i and depends exclusively on the atom type. We 
consider the simple case of a disordered crystal, made of 
one kind of atoms, characterized by the scattering fac- 
tor /, and containing impurities of scattering factor //. 
To understand the role of the scattering factors, let us 
start by evaluating the case of fixed atoms {ui = 0). We 
obtain: 



FIG. 2: Example of 1-dimensional crystal. The position of 
the ith-aiom is given by the equilibrium position Ri and the 
displacement Ui. In the figure the black circles are the host 
atoms with scattering factor / and the white circle represents 
an impurity with scattering factor // 



and the symmetry of such peaks allows to fully charac- 
terize the structural properties of the CDW. In this sec- 
tion, we isolate the different terms which contribute to 
the satellite peaks, and study their symmetry properties. 
This discussion is general and model independent. 

The expression for the total diffraction intensity of a 
crystal ist^ 

^('?) = ^ E e-*^^^'-^^-^ (/,/,e-«(— .)) . (14) 

As sketched in Fig|21 Ui is the atom displacement from 
the equilibrium position Ri = ia, with a indicating the 
lattice constant. (...) denotes the double average over 



J 



K 



where A/ — fi — f and / is the average scattering factor. 
The usual Bragg peaks, corresponding to the reciprocal 
lattice vectors K, arise from the first term in H15|l . while 
the second term is responsible for a constant background 
intensity, called Laue scattering, due to the disorder. 

Now we move back to the general case Mi ^ 0. We are 
interested in the behavior of the scattering intensity I{q) 
near a Bragg peak {q ^ K). Since \5q\ = \q — K\ <^ K , 
we can take the continuum limit i r and obtain from 



Hq) = 



-i(5(j(ri— r2) 



fri Ir- 



-iSq(u{ri) — u(r2)) 



(16) 

Here /^^ - yb^ / rf'^rid'^rs and /. = / + A/a'^/2E(r). 
Assuming that in the elastic approximation the displace- 
ments remain small (uj Ri), one can expand H16|) as a 
power series of Kuq. Developing up to the second order 
we geifi^ : 



I{q) ^ Id + L+ ^tripl, with 



(17) 



Id-- f q I e 



-^**("('"^)-"("=»(M(ri)M(r2)), h = -iqAf f e-**«("(''i)-"('^^» (jE(ri)M(r2) - S(ri)M(r2)) 



ipl 



-iqAfa" 



/j:(^ri)i:{r2){u{n) - ^(ra)) 



While the contribution 1^ represents the intensity due to 
the atomic displacements alone, the contributions /a and 
Aripi are generated by the coupling between the disorder 
and the displacements. In the following we consider only 
/a and Id', the term /tripi is evaluated in the appendix IXI 
where we show that it is smaller than the other two. 

In a pure system (E(r) = 0) we expect that only the 
first term can be different from zero. Referring to the case 
of a CDW without disorder and neglecting the thermal 
fluctuations, we can replace in (|17|) the displacement term 



(|13|l with (f){r) = const. This gives: 

/d(g) = /VE^(«±Q-^)- 

K 



(18) 



The presence of two satellites around each Bragg peak is 
one of the most clear experimental evidence of a CDW. 
In a pure system these satellites are symmetric and 
without any broadening. To interpret the experimental 
flndingsiSiiLiS, which reveal the presence of asymmetric 
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K-Q K K+Q 



FIG. 3: Sketch showing the symmetry properties of the two 
different terms contributing to the sateUite intensities (the 
Bragg peak at K is not shown). Top: term Id- The intensity 
is symmetric with respect to the Bragg peak. Bottom: term 
/a. The intensity is antisymmetric respect to the Bragg peak. 

peaks, we need to account for the effect of impurities. 
In particular, we want to explain not only the measured 
intensity asymmetry (lA) between the two satellites, but 
also the profile asymmetry (PA) of each peak, which is 
measured in strongly doped samples. 

The symmetry properties of the terms /a and 1^ can 
be determined by considering the lattice displacements 
and the disorder expressed in terms of their Fourier com- 
ponents: 

u{r) = [ e-''"'uq 

JBZ 

= / e-^'-S,. (19) 

JBZ 

It easy to obtain: 

h{q) = fq\us,u^sq) (20) 
h{q) = 2a'^/2qA/7(/m[E_5,UA-,]). (21) 

The two prefactors, / and Affq vary slowly in q 
and one can assume they are constant in the neighbor- 
hood of the reciprocal lattice vector K. Due to the fact 
that (|21|l is an imaginary part, we can deduce: 

hiK + Sq) = hiK-Sq) 

h{K + Sq) = -h{K-6q). (22) 

In Fig. 121 we sketch the following symmetry properties 
of different terms: /a generates two satellites symmetric 
with respect to the Bragg peak, while /a gives antisym- 
metric contributions. We conclude that the lA is due to 
/a while the PA is not excluded in both terms, in par- 
ticular we have a mirror symmetry for I^- Plugging the 



displacement form (|13|l given by the FLR model in l|20|l . 
we get 

h{K + Q + k) = ulfK^ [ Cd(r) , (23) 

J r 

where ~ J d'^f and Cd{r) is the positional correla- 
tion function 

Cd(r) = (^e*(<^(i)-0(-i))^ . (24) 

Examining this equation one is able to connect the sym- 
metry properties of the satellite peak profile with the 
symmetry properties of the system. In particular, if 
Cd(r) is an even function in r, the Id term cannot show 
a PA. On the other hand Id{q) is a real function so the 
latter condition is equivalent to require the presence of 
a (j) —(f) symmetry in the system. We will make use 
these observations to interpret the experimental findings 
in section Hvl 

III. THE REPLICA METHOD 

In this section we calculate the different terms of the 
development of Eq. l(T7|) using a Gaussian variational 
approacbiiiS^. We consider first the isotropic case, to 
include later the corrections due to the Coulombian in- 
teractions. The physical interpretation of these results is 
presented in more detail in section Hvl 

We consider a FLR Hamiltonian H = H^i + iJdis where 
the second term is (|ll|l and the first one is described ei- 
ther by l(2Jl, in case of a short range elasticity, or by © 
taking into account the effect of long-range Coulombian 
interactions. We first perform the average over the disor- 
der using the replica techniques. The replicated Hamil- 
tonian is: 

i?eff - ^ci - / d'r^ E ('^«(^) - ^"(r)) , (25) 

a '' a,b 

where T is the temperature and the sum over the n repli- 
cas has to be considered in the limit n ^ 0. We observe 
that the system is — > —4> invariant. This means that 
the FLR model with a Gaussian disorder cannot generate 
satellite peaks with a PA. Finally, we stress that, moving 
from the original Hamiltonian to its replicated version 
we also need to change the correlation functions contain- 
ing explicitly the disorder. In particular, using l(T^. the 
terms Id and /a obtained in H17() become 

Id = Fd / e-»*'('-i-'-^)V,,V,,(pi(ri)pi(r2))eff 
I.-fJ e-^^^(---)^Il^^(p„(ri)p.(r2)),, 

a,b 

where Fd = and = "'"°qtvo^^^ - Replacing 

the backward scattering term of the electron density and 
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performing an integration by parts we get to the form: 



Id = f q u 



f e-'^ir y^-^Qr + ^.c] Cd(r) , (26) 

J r 



J r 

a 

1 " 



c.c.i aw 



gZ(0a(§)-0a(-i)) 



off 



gi(0a(§)-06(-5)) 



off 



are the positional correlation functions controlling the be- 
havior of each contribution. To obtain this result we have 
applied the standard decomposition in center of mass R 
and relative r coordinates (ri = i? + § and 7-2 — R~ 
Since u varies slowly at the scale of the lattice spacing, 
we performed the integration over R. We notice that 
(|25|l reveals clearly the presence of two peaks situated 
ai q = Q + K and q = K — Q. In particular, as ex- 
pected, the contribution to the two satellites of the dis- 
placement term Id is symmetric, while the one of /a is 
antisymmetric. The sum of these two terms leads to the 
lA experimentally observed. 

Following the same method used to study the flux lines 
in presence of a weak disorder—, we can evaluate the dif- 
ferent terms in H26() . We look for the best trial Gaussian 
Hamiltonian Hq — JqG~^{q)(j)a{q)4>b{—q) in the replica 
space which approximates (|25|l . The G~^{q) is the n x n 
variational matrix. Without loss of generality, this ma- 
trix can be chosen of the form = cq^Sab — crab- The 
connected part is defined as = ^i,G~j^. By mini- 
mization of the variational free energy we derive that Gc 
is given by the bare elastic propagator. In the isotropic 
case we write 

= cq^ . (27) 
For a long range elasticity, in c?-dimension, it follows: 



G-'^c^^ + cq\ 
Finally, the parameters Gab are given by: 

D Bat,('-=") 



(Ta^b = j;e 



where 



2T 



BZ 



G{q) - Gab{q) COS qr 



(28) 



(29) 



(30) 



G is the diagonal element of Gab- In the Gaussian ap- 
proximation, the positional correlation functions become: 



Cd(r) = 

1 " 



(31) 
(32) 



where B is the diagonal element of Bab- In order to 
simplify the notation we write H32|) as Ca(r) — x{''')Gd{r)- 
Finally, using (jSOJ, we get: 



Xir) = ^ E 



-'Tj(G{q)-G^t{g))cosqr 



(33) 



a,b 



Two general classes of solutions exist for this problem: 
while the first class preserves the permutation symmetry 
of the replicas (RS), the second class (RSB) breaks the 
replica symmetry. It has been shown^^ that the stable 
solution for d > 2 corresponds to the RSB class, while 
the RS solution remains valid at short distance. In the 
following we will refer directly to the d — 3 case. 



A. Replica symmetric solution 

We discuss first the RS solution which gives the cor- 
rect evaluation of the correlation functions at a distance 
smaller then Ra- Within this Ansatz the matrix G and B 
are defined by two of their elements: the diagonal values 
G, B and the off-diagonal ones, B — Ba^b and G — Ga^b- 
The simple algebra of symmetric matrix yields, for n — Q, 



G = G c{l + G cO a^b) : 

G G = Gr 1 



(34) 
(35) 



where from (|29(l and 10} , we deduce that aa^b ~ 
D/Texp{-'(t)^/2). 

We focus first on the short range elasticity. Replac- 
ing the form of Gc given in H27|l we easily calculate the 
displacement 

~ At 

37r^ i?a 

The first term take into account the thermal fluctuations 
and saturates to (fy^ at a large distance. The second term, 
due to the disorder, grows with a power law and it is 
responsible for the exponential decay of the positional 
correlation functions. To characterize the spectrum it 
remains to calculate 

1 



T 
2^2 



1 



(37) 



We can conclude that in the RS scenario the two po- 
sitional correlation functions decay exponentially fast, 
moreover Ca has a power law extra factor. 

Before calculating the stable RSB solution, we evalu- 
ate the scaling behavior of the same objects in the case 
of an unscreened Coulombian potential. The only change 
consists in taking the connected propagator given in H28() 
to determine the physical quantities. It is instructive to 
discuss first the general d-dimensional case. The dis- 
placement takes the form 



B{r) 



a.b 



qdqx 



q 



2d-2 



:[l-cos(gr)]. (38) 
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A similar integral was discussed in Ref. '29\. A general 
remark is that scales as The scaling behavior 

of these integrals is determined by small q's, for this rea- 
son we can neglect the Qx dependence in q. The strong 
anisotropy along x and y can be studied performing the 
following substitution: 



V = xqx 
z - 



t = qy 
c X 

d+l ■ 



Cly 



We obtain 



B{r) 



+2d-2 



[(f) 



(39) 



:[1 -cos(t + i;)] 



Hi{z)., 



(40) 



where iJi(O) = const, and Hi{z ^ oo) cx 2(9-3<i)/(d+i), 
It easy to check that, for d = 3, there are only logarithmic 
divergences: 



B{x^O,y) 
B{x,y^O) 



D 



(27r)2Vc3^ 
D 



log(Ay) , 



log(A^a;) . 



This result is a clear evidence that, because of the long 
range interactions, the system is more rigid and the crit- 
ical upper dimension becomes d = 3, in contrast with 
the result d = 4 for the short range case. To confirm 
this statement we evaluate the Fukuyama-Lee length by 
imposing B(r) = (27r)2: 



Ra(x = 0,?/) 
Ra{x,y 0) 



A; 



(41) 



The exponential law is characteristic of the upper critical 
dimension and is an extrapolation of the power law. The 
equivalent form ljB7|l has been derived for an isotropic 
short-range elasticity in d = \S0^3i ^ j-j. remains to deter- 
mine Eq. (|37|l . At low temperature we write 



d'^qq 



d-l 



COS qr 



d+l 



5-3d 
5-3d 



ciqj + cq 
dv t 



Y cos(t 



/CCl 



(42) 



where iJ2(0) = const, and H2iz ^ oo) cx 2(5-3<i)/(d+i)^ 
In d = 3 a straightforward calculation confirms this scal- 
ing behavior: 



Xix,y = 0) 
X{x = Q,y) 



1 



IGncx 
1.1 



We observe that the anisotropic scaling a; ~ j/^ is always 
verified. Since for d = 3 the RS solution is unstable, to 
obtain the physics at large distance one has to consider 
the RSB method. 



B. Replica symmetric breaking solution 

Within this scheme, the off diagonal elements of Gab{q) 
are parameterized by G{q,v) where Q < v < 1. The 
saddle point equation becomes 



a{v) 



D 



(43) 



We look for a solution such that (t{v) is constant above 
a variational breakpoint Vf.. This can be done in a easy 
manner by recasting the equations in terms of a new vari- 
able 



[(t](u) = va{v) — / du(j{u) 



(44) 



It is not difficult to show that [o']'(?;) = va'{v). We refer 
to appendix^ where we summarize the previous results 
for [a] and we calculate its form in the case of an un- 
screened Coulombian elasticity. As a first step one uses 
the inversion rules of hierarchical matrices^S , Eqs. H34f) 
and H35|l become 



G = Ge(l 

G-G = 



dv 



r) 



Gc' + [aY 
1 r .. I'it) 



G7' 



dt- 



{G^' + [a]{t)Y\ 



(45) 
(46) 



where E = [cr](wc) is a variational parameter, whose ex- 
pression is determined in appendix IbI 

Starting from the short range case we calculate the dis- 
placement using Eq. (|45|l . replacing the expression ljB3| 
for [a]ii: 



Bir) 



hi 



47r' 

IBZ 

2 log(Ar) . 



[1 — cos qr] 



(47) 



The logarithmic behavioB^iiii of 1)47(1 is controlled by small 
V {v < Vc)- Values of v above the breaking point {v > Vc) 
give the small distance contribution. To fully character- 
ize the spectrum it still remains to evaluate xl*") in. the 
RSB scenario. 



Giq)-G{q,v) 



2 

d2 



dt- 



(48) 



where the parameters I and Vc are given by (|B4|1 . By 
integrating (|48|l over q and with some manipulations, we 
get: 



X[r) = Y 



1 - 



dz cxp — Btt 



t 



(49) 
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FIG. 4: Ratio between the RSB and RS solutions of x{^)- 
At large distance this ratio tends towards a constant value 
b, with b ~ 0.018. This means that the RSB solution affects 
only by a multiplicative factor. 



The low temperature behavior {I ^ Ra) of this term is 
sketched in Fig. 0] 

Finally we study the RSB solution for the Coulombian 
elasticity. Replacing ljB8|) in (|45l) , the diagonal correlator 
becomes 



G 




(50) 



The latter form is characteristic of the upper critical di- 
mension. Inserting Eq. (|50|l in Eq. (|30|l , taking ci = c — 
1, and employing the usual substitution H39() we get: 



B{r) = S2 dt 



dv 



r[l -cos(^ + ^;)] 



This equation leads to the same conclusions discussed for 
the isotropic upper critical dimension. The asymptotic 
displacement is thus given by 



Bix) ~ log(log(A,x)) . 



(51) 



As in the case of the short-range elasticity, it can be 
shown that the RSB solution does not affect the asymp- 
totic power low behavior of x(^)- 



IV. PHYSICAL DISCUSSION 

In this section we summarize the results obtained in 
the previous sections and compare them with the exper- 
imental findings. 

We start by evaluating for the short-range elasticity 
model the positional correlation function Cd(r), defined 
in Eq. (|24|l . This function is the analog of the correla- 
tion function determined for vortex line systemgii. It is 



FIG. 5: Behavior of Cd{r). The thermal contribution satu- 
rates to exp(— 0y/2), the RSB solution is valid beyond Ra, 
while the RS one is valid below Ra- 



well known that for > 2 it exists a finite temperature 
(the critical temperature Tc) below which thermal fluc- 
tuations are prevented from disordering the system. In 
fact, the direct calculation of the positional correlation 
function in d = 3 yields: 



C. 



therm. 



(r) 



e 2 



.(1- 



(52) 



with Si(Ar) — J^^ dtsint/t. As shown in Fig. |S1 the 
correlation function saturates after a few lattice parame- 
ters to a non-zero value, which witnesses the presence of 
long-range order in the system. 

However, the quenched noise originated in the impuri- 
ties is still able to destroys the long range order, even in 
d — 3 (see Fig. (S)). In this circumstances, the more tra- 
ditional scheme^^ describes the system as organized in 
ordered domains (called 'Fukuyama-Lee-Rice domains'), 
characterized by the average size Ra- Beyond this charac- 
teristic length, the CDW dislocations become dominant 
and any order disappears. This scenario is captured by 
the RS solution. From ^ we obtain: 



(53) 



As a direct consequence, we expect to find Lorentzian 
satellites, whose half-width at half-height is of the order 
of R~^. However, the RS approach is unstable at large 
distances and the correct solution is given asymptotically 
by the RSB. Within this assumption, the correlations 
decrease following a power law: 



CfSB(r) = C'f'=™-(r)(— )". 



(54) 



Applying the variational approach , one finds, from 
(07)), rj = I. The corresponding quasi-ordered phasei2*ii, 
called Bragg glass, is characterized by an infinite cor- 
relation length and the characteristic size Ra represents 
now the crossover between the RS and RSB solutions. 
One can give a simple physical interpretation of the two 
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FIG. 6; Intensities of the different contributions to satellite 
peaks. The more divergent term, Id, is symmetric. 7a is 
antisymmetric. In this figure we consider a repulsive potential 
and A/ > 0. 



identified regimes, by observing that at the scale Ra the 
phase distortions are of the order of the CDW period 2tt. 
This means that for distances smaller than Ra the devel- 
opment of the Hamiltonian H25|) is allowed and leads to 
the RS solution. For distances larger than Ra, instead, 
the phase feels its periodic nature and this trivial devel- 
opment of the Hamiltonian is no more valid. 

In order to determine the intensity and the shape of 
the two satellites we need to evaluate all the terms con- 
tributing to the development (|17|l . In particular, we con- 
sider the displacement term Id, and the asymmetric term 
/a, arising from the coupling between disorder and dis- 
placement. In the literature this latter term was pre- 
viously estimated by means of modelsi^iiSiiS which de- 
scribe the pinning imposing a constant value (jjQ on the 
phase in proximity of each impurity. According to these 
approaches, the observed satellite asymmetry is a clear 
signature of strong disorder—. Thanks to our more ac- 
curate calculation, we found that the term /a is non-zero 
also in case of weak disorder and it gives rise to a diver- 
gent contribution similar to the one stemming from I^. 
Using (H?)), (|2ZIl and our final results read: 

I^{K + Q + k) = J^KW f e-'^^^i^Y, (55) 

J r ^ 

UK + Q + k) = -2^^jAf l^KuJ e-^^^{^r^-^. 

y Ra Jr r r 

After computing the d-dimensional Fourier Trans- 
forms, we conclude that both terms are divergent: in 



FIG. 7: Dotted line: divergent peak Id{<l) ~ for a perfect 
experimental resolution. Full line: peak convolved with a 
finite experimental resolution according to 15611 . The height 
of the peak is ~ ^^Ra, the half-width half-height ~ is 
also indicated. Dashed line: an approximate form for the 
peak = r/''d{0){l ~ q^e/i). 



would be dominant, as Cd(r) decays to zero less rapidly 
than Ca(r). However, in a real measurement the diver- 
gence in is always cut by the finite resolution and 
both terms have to be accounted for, as the prefactor 
Kuq of the less divergent term /a is larger. 

At this stage, following the analysis of the experimental 
data of Ref. pjj concerning neutron diffraction spectra, 
it is interesting to discuss the role played by the experi- 
mental resolution in determining the peak shape. With 
the low resolution achieved by means of neutrons to de- 
scribe a vortex lattice, in the above cited experimenfei^, 
it is possible to determinate the intensity spectrum only 
along one direction, after performing an integration on 
the other two directions. Instead, the much higher reso- 
lution reachable in X-rays experiments, is in principle ad- 
equate to perform the whole 3-dimensional Fourier trans- 
form of the spectrum. For concreteness sake, we assume 
a Gaussian resolution with variance where ^ > Ra 
(the opposite case is not interesting as the resulting peak 
shape is affected only by the resolution). We consider 
first the behavior of the direct term Id in case of a RS 
solution. The peak obtained from the C^^{r) function 
is essentially independent of the resolution and has the 
shape of as a squared Lorentzian of height cx R^ and 
half- width cx 1/Ra- The profile drastically change if we 
consider instead the correct RSB solution. Setting ?/ = 1 
we can find an analytical expression for the experimental 
peak: 



particular, /d oc 



and /a oc 



This result, sum- 



<7° 

marized in Fig. |2| is a clear sign of the quasi-long range 
positional ordered phase. In this particular case, the peak 
a.t K + Q is smaller than the K — Q one, since the specific 
interaction between the impurity and the CDW is repul- 
sive (we would have the opposite asymmetry in case of 
an attractive interaction). We observe that for an ideal 
experiment with infinite resolution the symmetric term 



k) ^ Rl 
= (27r)3/2 



(56) 



k 



k^ 



-er/.(^). 



where erfi(x) = —ierf{ix) is the imaginary part of the 
error function. This expression gives a non-divergent 



—2 _ 



peak, shown in Fig. [T] whose height is Air K'^uq^ R^^'^ 
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FIG. 8: Sketch of the experimental finding*—. The two satel- 
lites present an I A in agreement with our predictions, but also 
an evident PA. 



and its half-width ~ 1.5034/^. We remark that the 
peak decays with the characteristic law only for 
< g <C while a wide region around the max- 
imum height is dominated by the effects related to the 
finite resolution. In Fig. [7| we also report the function 
= /r^(0)(l - It is possible that, in 

a real experiment, the signal-noise ratio is low enough to 
hide the behavior. In this specific case, to put in ev- 
idence the existence of a Bragg glass regime one should 
vary the Fukuyama-Lee length -Ra^*"- In fact, if the peak 
shape is Lorentzian, varying Ra produces a change both 
in the width and the height of the peak, whereas if the 
peak follows a power law, only its height is changed due 
to a variation of Ra, while the width is fixed by the res- 
olution. 

In absence of screening, the long range Coulomb in- 
teractions become important: the system is more rigid 
and the upper critical dimension is shifted from d = A 
to d = 3. As a consequence, the satellite peaks become 
more divergent. In particular, we find that the symmet- 
ric term Id goes asymptotically as instead of q~^ . 
Moreover, the strong anisotropy between the longitudi- 
nal direction (x) and the transversal ones (y) leads to 
an anisotropic scaling of the correlation functions. We 
verified that in d = 3 the dependence on x and y of the 
correlation functions (as clearly shown in 1021) a-nd of 
the characteristic length scale (|41ll respects the relation 
/(cc, y = Q) ^ f{x = 0, y'^). As a result, in this regime, we 
expect more divergent and more anisotropic peak shapes 
in comparison to the ones observed in the short-range 
case. 

On the experimental side few detailed diffraction spec- 
tra are available at the moment. An example of proto- 
type systems are doped blue bronzes?^^, where the dis- 
order is introduced by partially replacing Mo^+ for the 
non-isoelectric V^^. In this kind of system, the interac- 
tions between the CDW and the impurities are repulsive 
and A/ is negative. The shape of the observed spectrum 
is sketched in Fig.|Sl The sign of I A is consistent with our 



prediction as well as with the predictions of the strong 
pinning approach. However, we also remark the presence 
of PA of each single satellite. The mirror symmetry of 
the shape of the two peaks suggests that the PA is origi- 
nated in the correlation function Cd{r). As we have pre- 
viously discussed, the correlation function obtained from 
the FLR model with Gaussian disorder is real and pro- 
duces two peaks with profile symmetry. One can wonder 
what is the effect of non-Gaussian disorder. To this pur- 
pose, we studied a 1-dimensional model with a binomial 
distribution of impurities. If we restrict to the forward 
scattering term in the development (Q, the problem is 
exactly soluble. We find two Lorentzian peaks without 
PA, centered around a vector KzLQ shifted with respect 
to K±Q as a consequence of the presence of odd moments 
in the disorder distribution. It remains to investigate the 
role possibly played by the amplitude fluctuations. The 
short correlation lengths extracted from the experimen- 
tal data suggest that this particular effect is most likely 
to happen in the strong pinning regime, whereas our cal- 
culations concern the weak pinning limit. The authors 
of Ref. (d) indeed justify the PA with the presence of 
Friedel oscillations and hence with the presence strong 
fluctuations of the amplitude of the condensate, at least 
in the neighborhood of the impurities. 



V. CONCLUSIONS 

Summarizing, we determined the shape and the inten- 
sity of the satellite peaks characterizing the spectrum 
of a pinned charge density wave. We analyzed in de- 
tail the case of a weak and collective disorder, when the 
Fukuyama-Lee-Rice model is justified. We considered 
both the short-range elasticity as well as the long-range 
elasticity generated by an unscreened Coulomb interac- 
tion. In both these cases, we found divergent peaks dis- 
playing intensity asymmetry. The divergent nature of 
the peaks is, as it was discussed, the clearest sign of a 
Bragg glass phase. Moreover, the long-range elasticity, 
when present, is responsible for a larger anisotropy and 
a stronger divergence. Let us stress that the calculated 
sign of the intensity asymmetry is in agreement with the 
experimental data. We discussed the role played by the fi- 
nite resolution of the experimental setup, calculating the 
convolved shape of the peaks, where the divergence is cut. 
From these observations, we illustrated possible methods 
to reveal experimentally the presence of the Bragg glass 
phase. Concerning the asymmetry of the peak profile, 
we showed that, on general symmetry grounds, it is not 
expected in the weak pinning regime. We conjectured 
that its observation in a recent experiment is likely due 
to the strong pinning present in the measured system. 
Finally, we observe that the profile asymmetry may hide 
the power law behavior of the satellite peaks. It would 
thus be highly desirable to dispose of measures in less 
disordered systems where one can expect a Bragg glass 
behavior, e.g. using isoelectric impurities. 
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APPENDIX A: THE TRIPLET TERM 

In this appendix we discuss the behavior of the term 
^tripi of the intensity development (|17|l . This term was 
only conjectured to be negligibl o^^i?"^ , before tha^ a di- 
rect calculation of it was performed. 

We start analyzing the symmetry properties of /tripi 
showing it has the same symmetry of /a. Next we eval- 
uate /tripi within the RS scheme. Independently of the 
elasticity range, we find that the difference between the 
two contributions is given by a simple g-independent mul- 
tiplicative factor. 

Applying a Fourier Transform (|19|l we obtain 



/tripi(i^ + Sq) - ~2tKAfa'' [ (S^-,!], 

JBZ 



(Al) 

This equation gives a non-zero contribution only if we 
consider higher harmonic terms of the electron density 
(HJ: p{x) ~ pi\^p\cos{Q{x + (t){x))) + pi\ip\cos{2Q{x + 
(j){x))). The two satellites take the form: 



Itrip\iK + Sq) = -2iK A f a {'Esq[^-2SqUSq + '^SqU-2Sq]) 

At this stage it is evident that /tripi has the same symme- 
try of /a. Performing the integration over the Gaussian 
disorder by means of the standard replica techniques we 
obtain 

/tnpi - AfquoNja^D I e'^''^^ [e-'^'' - c.c] CtHpi(r) 

J r 

Using the usual decomposition Ctripi(7') = "xtripK'^) 
we write 

A_Y. [e-^^-^^I'^-c., (A2) 

„-Tjgj,[(G-G„i,) + (G6<:-G„6)]cos<7r 



Xtripi(r-) 



We introduce the replica symmetric Ansatz. It is easy to 
check that for n = 

~ ^ ^ ^gfcc — Aaaa — '^^^ (^Aaai}-\-Aai,a-\-Ai,aa)-'r2 ^ ^ Aabc 



abc 



a^b^c 



Using this relation and (|35|l we can evaluate Xtripi(''). 
In order to simplify the notation we recall that (/i^^^ — 
2T Jr 7;Gc and x{r) = [1 - exp{-T J^^ GcCOsqr)]/T. 
From (IA2I) we obtain 



Xtripl{r) 



X{r) 
T 



1 - ( 2sinh(T [ Gccosqr) + 1 

\ Jbz 




K-Q 



K+Q 



FIG. 9: Intensities of the antisymmetric contributions to 
satellite peaks. We remark that /a S> /tripi. 



Because we are interested in the long distance behavior 
we remark that J^^ GcCosqr — > 0, whenever r — > oo. 
Developing up to the first order we get 



Xtripiir) x(r) 
- Xir) 



1 - 



T 



2T 



(A3) 



where the last step is valid at low temperature. In fact, 
in this regime Ctripi(?') and Ca(r) have the same behavior 
in r. If we take a spherical cut-off A = 2TT/a, where a 
is the lattice space, the integral (j>'^/{2T) — (cTra)~^ is 
independent of the range of the elasticity (for a ^ 0). 
We can compare the two terms /tripi and /a, by writing 
the ratio between the two intensities 



'tripi 



A/ 



Nia 

Ra 



(A4) 



For weak disorder Ra ^ a; it follows that /a » /tripi- 

The evaluation of this term in a more accurate RSB 
approach is very complicate because Xtripi {f) involves the 
sum over three replicas. However, in analogy with %('') 
we can argue that the RSB solution does not affect the 
asymptotic power low behavior of Xtripi (t) ■ In Fig. we 
summarize our result taking the correct RSB behavior 
for B. 

We conclude that the triplet term renormalizes the 
prefactor of /a without changing the power law behavior. 
In particular, when A/ > 0, /tripi enhances the asymme- 
try between the two satellites, while when A/ < 0, /tripi 
slightly decreases the antisymmetric contribution. 



APPENDIX B: CALCULATION OF [a] 

In this appendix we determine the variational function 
[(t]{v) defined in Eq. (|44|l . We start form the saddle point 



12 



equation H43|) . From H30() we know that 



B{r = 0,v) 



G{q)~Giq,v) 



(Bl) 



BZ 



The integral in the momentum space is performed in the 
Brillouin zone. To simplify the analytical form of our in- 
tegrals the ultraviolet cutoff, A, is taken equal to infinity 
whenever it is possible (i.e. whenever the integrals are 
ultraviolet convergent). Inserting H46() in the previous 
equation and taking the derivative of H43|l leads to the 
equation determining [a]: 



df^q T 



= 1 



(B2) 



Solving the integral for 2 < d < 4 (A — > oo) and deriving 
again, one gets 



[a]{v) T.{v/vcf/'^ forw < Vc 
~ S for V > Vc ^ 



(B3) 



where 9 = 2 — d. The values of the breakpoint Vc and 
S = cl~'^ determine the crossover between a short dis- 
tance regime {x I), where [a] is constant and the RS 
solution valid, and the asymptotic regime (a; > I), where 
the physics is determined by the small v behavior of [a]. 
Using Eqs. (|B2p and H43() . after some manipulation it is 
found for d = 3 



I = 



T 

8lc' 



(B4) 



We observe that the crossover between the two regimes 
is I ^ Ra in agreement with the dimensional result (|12|l . 

For d — 4:, and more in general at the upper critical 
dimension, the integral in l|B2|l has a logarithmic ultra- 
violet divergence. As discussed in Ref. [S^l, the behavior 
of [a] when v is small is not described by a pure power 
law. Starting from (|B2|) for d = 4 we get 



a{v) 



d'^q 



T 



(27r)4 [cq^ + H]2 

Jo [q^ + 1]' 

S4<j{v)T cA\ 



(B5) 



Where Sd = 2'^-'^n-'^/'^ /T{d/2) is the angular integration 
in d dimension. Defining A = 2c^ /{S4T), we obtain after 
one more derivative 



Av 



Av 



(B6) 
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FIG. 10: Continuum line: integral over q in <B8ll as a function 
of Qx. Dotted line: limiting behavior I[qx] Tr/Sqx- Finally 
I[0] = 0.5 



This result is valid up to log- log corrections. From 
Eqs. ljB5(l and H43|l we can estimate the crossover length 
I. In the short distance regime, where [a] = cl~^, it turns 
out: 



A-ie^ 



(B7) 



At this stage we can discuss the case of the Coulombian 
long range elasticity. Starting from Eq. ljB2|) and using 
Eq. ||2H1) for Gc, we write 



d^q 



T 



f^" dqx 



To solve this equation we consider the physical case, 
where A^: Q, A ~ 27r/a and a is the lattice space. 
In this limit A^, ^ A so we can assume A — > 00 and solve 
the integral 



T(7{v) 



"'0 



47r^Vc?ci 7o 
q^dq 



dqxl[qx 



{ql + q^ + iy 



The behavior of I\qx\ is sketched in Fig. EH In conclu- 
sion, we obtain 



AlrV 



AlrV 



log 



2 Awv 



(B8) 



where Air — IQir^/c^Ci/T . This equation is equivalent 
to the one found for an isotropic system at the upper 
critical dimension. 
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